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Abstract 

Microscopic robots could perform tasks with high spatial precision, such as acting on precisely- 
targeted cells in biological tissues. Some tasks may benefit from robots that change shape, such 
as elongating to improve chemical gradient sensing or contracting to squeeze through narrow 
channels. This paper evaluates the energy dissipation for shape-changing (i.e., metamorphic) 
robots whose size is comparable to bacteria. Unlike larger robots, surface forces dominate the 
dissipation. Theoretical estimates indicate that the power likely to be available to the robots, as 
determined by previous studies, is sufficient to change shape fairly rapidly even in highly-viscous 
biological fluids. Achieving this performance will require significant improvements in manufac¬ 
turing and material properties compared to current micromachines. Furthermore, optimally 
varying the speed of shape change only slightly reduces energy use compared to uniform speed, 
thereby simplifying robot controllers. 


1 Introduction 


Metamorphic robots [l][6||l6, 35}[36][43j can change their shape to operate effectively in variable un¬ 
structured environments, e.g., for searching collapsed buildings. Extending this concept, ongoing 
progress in reducing robot sizes toward that of biological cells 38 raises the possibility of micro¬ 
scopic metamorphic robots. Such robots could be particularly useful in medical applications 13 
by operating in biological fluids. Examples include robots matching the shape of cells to provide 
temporary scaffolds during tissue repair, and extending probes into spaces, e.g., between or into 
cells, that are too small or fragile for the whole robot to enter. 

Microscopic organisms and robots face significantly different operating constraints than larger 
ones 32 . In particular, their high surface to volume ratio means surface forces, such as viscous 
drag and internal friction, dominate the behavior of microscopic robots. Surface forces are espe¬ 
cially significant for shape change, which involves external surfaces moving through the surrounding 
fluid and internal surfaces moving past each other to actuate the shape change. Moreover, the 
need to scavenge power from their environment, rather than relying on tethered power, limits the 
power available to the robots. Thus, in addition to the challenges of manufacturing and operating 
microscopic robots in general, power could constrain how rapidly metamorphic robots change shape. 

This paper evaluates power required for microscopic metamorphic robots to change shape in bio¬ 
logical fluids, including the effects of both external viscous drag and internal actuator friction. This 
contrasts with one application locomotion - where periodic shape changes have been extensively 


studied 29 by focussing almost exclusively on external viscous drag. Specifically, the remainder of 


this paper discusses power dissipation for microscopic metamorphic robots, and then evaluates the 
dissipation for three prototypical examples: a robot expanding its length, a robot extending a thin 
probe and robots coming together to form an aggregate structure. 
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scenario LOW high 

density p 1000kg/m a 1000 kg/m 3 

viscosity 77 10 ~ 3 Pa-s 10 Pa • s 

kinematic viscosity v = 77 /p 10 _ 6 m 2 /s 10 _ 2 m 2 /s 

Table 1: Parameters for biological fluids similar to water and cell cytoplasm. 


2 Power Dissipation During Shape Change 


In biological tissues, viscosity varies by orders of magnitude while density is roughly the same. Thus 
as example scenarios we consider a wide range of viscosity but constant density. These scenarios are 
changing shape in 1 ) a low viscosity fluid comparable to water and blood, and 2 ) a high viscosity 
fluid 10 4 times more viscous than water, which is typical for mucus or cell cytoplasm 13 . Table [l] 
gives the fluid parameters. 


Power use depends on the robot’s geometry, the speed with which it changes shape, the viscosity 
of the fluid it operates in and the friction of its internal mechanisms. To isolate the effect of changing 
geometry, we define characteristic robot size (do) and speed of its shape change (uo). We take these 
values to be fixed for a particular scenario, thereby separating effects of shape from those of overall 
size and speed. 


The definitions of these characteristic quantities are somewhat arbitrary. For example, the size 
of a spherical robot could be its radius or diameter. Moreover, a robot may change size and have 
parts moving at different speeds as it reconfigures. Thus vq could be the maximum or average speed 
of all moving robot parts. With a specific choice of these characteristic quantities, the analysis 
presented here isolates the effects of shape by defining dimensionless geometric factors relating 
power dissipation to shape. The dependence of the results on do and vq indicate how power varies 
with respect to overall size and speed. 


The following analysis uses the characteristic size and speed to determine the suitability of 
simplifying approximations. This is suitable provided all sizes and speeds during the reconfiguration 
remain close enough to these characteristic values for the approximations to remain valid. Thus, 
for instance, we do not consider robots that change shape by orders of magnitude in size or have 
similarly large variations in the maximum speed of moving parts during the change. 


2.1 External Power Dissipation 


To determine power dissipated in the fluid surrounding the robot as it changes shape, we make two 
approximations to the fluid flow. First, for small objects moving relatively slowly in fluid, viscous 
forces are significantly larger than inertial forces. This leads to smooth, laminar motion called Stokes 
flow 


17 22 41 . This approximation is accurate when the flow’s Reynolds number 


Re = 


v 0 d 0 


(1) 


is small, where v is the fluid’s kinematic viscosity. 

The second approximation is the robot changes shape slowly enough that the fluid motion at each 
time is close to the static flow associated with the instantaneous geometry 
approximation requires a small value of the Womersley number 


22 . This quasi-static 


W = 



( 2 ) 


where T is the time over which the geometry changes. 

Since characteristic size and speed are somewhat arbitrary, so are the values of Re and W. 
Thus their order of magnitude rather than precise values provide a guide to the validity of the fluid 
approximations used here. 
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These two approximations give a simple scaling of energy dissipation with size and speed, allowing 
a numerical solution with one choice of parameters to also apply to other cases with small Reynolds 
and Womersley numbers. Specifically, the power dissipated has the form 

-Pfluid = ^fluid^doUg (3) 

where 77 is the fluid viscosity and hfi u [o is a dimensionless geometric coefficient depending on the 
robot’s shape and how rapidly it changes compared to v$. Thus, for given choices of do and vq, /ifi u id 
determines how power dissipation varies as the robot changes shape. 


2.2 Internal Power Dissipation 


Internal dissipation arises from the robot’s actuator mechanisms. Friction in the actuators is the 
major source of internal dissipation directly related to a robot’s shape change. Friction is a com¬ 
plicated property, depending on surface microstructure, lubrication, applied forces and operating 
speeds 25 39 . Robots made of precisely structured, stiff materials could, in theory, experience 


much less friction than current micromachines. Specifically, phonon scattering dominates the fric¬ 
tion of smooth surfaces separated by atomic-scale distances moving past each other at speed v well 
below that of sound in the material 91. This leads to power dissipation of the form 


Pfrict ion — ^-friction 


(4) 


where S is the area of the moving surfaces and Aifriction is a constant depending on the materials 
and their spacing. Thus, unlike the complexity of friction in micromachines with poorly defined 
geometry at nanometer scales, friction in smooth nanoscale surfaces has a simple dependence on 
area and speed. In particular, friction depends on the contact area S, in contrast to macroscopic 
systems with friction independent of the area [8j. 


For stiff materials, theoretical estimates of ^friction give values somewhat less than 10 3 kg/(m 2 s) 9, 
13 . We use this upper bound as a conservative estimate of internal friction for hypothetical actuators 


with stiff, atomically-flat sliding surfaces, i.e., we take 


fcfriction = 1000 kg/(m 2 s) 


(5) 


For comparison, the viscous friction for flat surfaces moving with relative speed v while sepa¬ 
rated by a layer of fluid with viscosity 1 7 and thickness d is r/Sv/d. This dissipates power ( r)/d)Sv 2 , 
so fcfriction = 77 /d. Extrapolating to distances of a few tenths of a nanometer gives fcf r i c tion ~ 
10 7 kg/(m 2 s) for fluids with viscosity similar to water. A similar value for fcf r i c tion comes from 
measured drag of micromachine rotors 7 . More precisely defined micron-scale surfaces can have 
lower friction 42 , but still significantly larger than the value of fcf r i c tion of Eq. |5j Thus the inter¬ 
nal dissipation we consider, while achievable in principle, is considerably smaller that of current 
micromachines. 


For the small robots considered here, sliding surfaces are a relatively simple structure. As shown 
below, the internal dissipation from sliding is relatively small compared to dissipation due to viscosity 
in the surrounding fluid, provided smooth materials with low values of fcf r i c tion become available. 
On the other hand, if improved material fabrication can only reduce friction to values an order of 
magnitude or so larger than that of Eq. [5j alternate designs might avoid having internal dissipation 
dominate the total power use. Examples of such designs include rolling bearings between surfaces 
or lubricants found in larger-size robots. Another possibility is using flexure joints, which do not 
involve sliding friction. 

We consider the situation where sliding friction dominates internal power dissipation. Other 
power uses, such as for the control circuits or losses in distributing power to the actuators, are not 
considered here since these forms of dissipation are not directly related to the changing shape. Thus 
we take / ’internal ~ Oriction ■ Wi t il Eq. i internal dissipation then has the form 


^internal interna I /'friction dgUg 


( 6 ) 
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where d q gives a characteristic area and /internal is a dimensionless geometric coefficient for internal 
dissipation. 


2.3 Energy Dissipated to Change Shape 

The energy used to change shape is in time T is 


E = 



(7) 


where 

P(t) = Pfluid(t) + Pinternal(i) ( 8 ) 

is the power dissipated at time t. The power varies with the robot’s shape through the geometric 
factors /ifluid and /internal in Eq. i and Eq. [gJ respectively. 

As an example, suppose the shape change involves a single degree of freedom, specified by a 
parameter f(t) monotonically increasing from f = 0att = 0tof = Fat time T. Uniform motion 
in terms of this parameter corresponds to f(t) = ( F/T)t. 

Both external and internal power are proportional to the speed squared, so the geometric factors 
/ifluid and /imternai are proportional to / 2 , where / = df/dt. Thus Eq. [ 8 ] has the form P(t) = fc (/)/ 2 
where k(f) is a proportionality factor depending on the robot’s configuration (through the parameter 
/), fluid viscosity 17 , and fcf ric tion- 

One application of this formulation is determining how to vary the speed of the shape change 
to minimize the energy dissipated. This amounts to Ending the function f(t) that minimizes Eq. [T] 
subject to the conditions /(0) = 0 and f(T) = F. The Euler-Lagrange equation determines the 
minimizing choice of f(t) as the one giving constant power use throughout the shape change. Thus 
the optimal rate to change the parameter is / oc l/^/k(f). That is, the robot should change more 
slowly in configurations with larger energy dissipation in inverse proportion to the square root of 
that variation. 

In the examples of shape change discussed below, dissipation has only a modest dependence on 
the changing robot shape. Thus energy-minimizing change has relatively little variation in / so that 
uniform motion is close to optimal. This is a useful practical observation: in these cases, there is 
little benefit from an elaborate control program to minimize energy instead of moving actuators 
uniformly. Since microscopic robots have limited computational capability, this observation allows 
using simple control algorithms without sacrificing much energy efficiency. 


3 Linearly Expanding Robot 


Example applications for microscopic metamorphic robots are 1) using a compact shape to travel 
through the bloodstream and expanding at an injury location to aid repair [14] , 2) temporarily 
increasing length to improve detection of chemical gradients [ll , 3) changing shape to alter effects 
of fluid drag and Brownian motion 12 , and, 4) forming close contacts with cells or other robots. 


This section examines an illustrative example of these applications: a robot that changes its length. 


3.1 Geometry 

Consider a metamorphic robot consisting of telescoping segments, allowing it to expand in length, 
as shown in Fig. [l] The robot consists of an outer cylinder, of length 2 L and radius L. The robot 
expands from both ends of the outer cylinder, with each side having n internal segments, each of 
length L. The innermost segment has radius r = L — ns and the rest are shells of thickness s. 

We suppose all segments extend simultaneously, with all segments continuing to move until the 
robot reaches its final configuration. Thus a single parameter specifies the robot shape, namely the 
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2 L 



2L(fn + 1) 


(b) 

Figure 1: Expanding robot with n = 5 internal segments on each side, with each internal segment 
extended a fraction / of its full length, (a) Telescoping cylinders with / = 0.5. (b) Cross section 
shown with thicker shells than parameters used for evaluating behavior (given in Table [2]). 
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expanding robot 



Figure 2: Coefficients for viscous drag, /ifluid, and internal dissipation, h internal) vs. extent of an 
expanding robot with n = 5 segments on each side and ratio of innermost to outer radius r/L = 0.75. 


fraction of extension for each segment, /. With uniform expansion from / = 0 to / = F in time T, 
/ = F/T and convenient choices for characteristic size and speed are the segment length and speed 
of the ends of the expanding robot: 


do = L (9) 

v 0 = nLf (10) 


3.2 Behavior 


To evaluate fluid drag we simplify the geometry to truncated cones extending from each side of a 
cylinder of radius L and length 2 L. The end of the cones move at speed nLf, and the speed along 
the side interpolates linearly from 0 at the base to nLf at its end. This simplification corresponds 
to the robot having a thin elastic outer layer covering the actuating shells. Using this geometry, 
we solve the steady-state Stokes flow with specific choices of size and speed to determine power 
dissipation as a function of extension. Dividing result by do^o gives Tifiuid, shown in Fig. [2j which 
does not depend on the size and speed used to solve the fluid flow. 

Evaluating internal power dissipation requires the sliding surface area and speed. The sliding 
surface area between segments i and i — 1 is Si = 27r(r + (n — i)s)L( 1 — /). The surfaces slide past 
each other with relative speed Vi — Vi-\ = Lf. Internal power use arises from the sliding friction 
between adjacent shells (Eq. [4]): 


Pi 


internal 


= 2 


i =2 


^friction c Lfr 


(ii) 


with the factor of 2 from the segments expanding in the opposite direction. With Eq. [6j this gives 

n — 1 /, , ,r\ , 

( 12 ) 


Tl — J. / V \ 

Paternal = 27t(1 - /) ^3 [i n ~ 2) + (n + 2)-J 


Thus ^internal decreases as the robot expands, i.e., as / increases, since each segment is in contact 
with less of its neighboring segments. Fig. [2] shows the behavior of p nterna i- 

With a fixed ratio between radius of innermost and outermost segments and fixed expansion 
speed of the ends vq, the dissipation from sliding friction decreases with n. I.e., using a larger 
number of thinner segments requires less power than a small number of thicker segments. Thus, 
the increase in sliding areas due to having more segments is more than offset by the slower relative 
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scenario 


LOW 

HIGH 

time to extend 

T 

1 ms 

Is 

speed of extending tip 

Vo 

2.5 mm/s 

2.5 gm/s 

Reynolds number 

Re 

2.5 x 1CT 3 

2.5 x 10" 10 

Womersley number 

W 

0.03 

10" 5 

fluid power factor 

VdoVo 

6.3 x 10~ 3 pW 

6.3 x 10~ 5 pW 

internal power factor 

^friction ^ 0^0 

6.3 x 10~ 3 pW 

6.3 x 10~ 9 pW 

energy used 

E 

1.5 x 10" 16 J 

1.4 x 10 -15 J 

glucose molecules 


70 

700 


Table 2: Expanding robots for the scenarios of Table |T[ with do = L = 1/rm, r = 0.75 ^m, n = 5 
shells, and fcf r i c tion = 10 3 kg/(m 2 s) for a robot expanding at a uniform speed from / = 0 to / = 
F = 0.5. The energy listed is for expanding at a uniform rate, given both in joules and the number 
of glucose molecule oxidations required to produce that energy, at 50% efficiency 20 . In both 
scenarios, fluid viscosity accounts for over 90% of the energy used. 


speeds between neighboring segments. This discussion does not include losses in the actuators: more 
segments require more actuators; and there is a minimum feasible thickness to the shells, e.g., to 
maintain required stiffness and have room for any required transmission of signals or power through 
the segment walls. 


3.3 Example 


Table [2] illustrates the behavior for the two scenarios of Table [lj The Reynolds and Womersley num¬ 
bers are small in both cases, indicating that quasi-static Stokes flow is an adequate approximation 
for the fluid behavior. 


Combined with the geometric coefficients shown in Fig. [2j these values indicate viscous drag is 
the dominant contribution to overall power use, and the power required for extending the robot in 
the times indicated is well below 1 pW. This is well within the power likely available to robots of 
this size, e.g., from chemical energy 13 20 . 

The geometric coefficients, shown in Fig. [2] have relatively little variation as the robot expands. 
This means the uniform expansion rate used in this example is close to the expansion that minimizes 
the energy. Specifically, the optimal expansion only reduces energy use by about 0.1%. 


4 Telescoping Probe 


A major capability of small robots is their ability to access small volumes, e.g., in biological tissue. 
Nevertheless, some locations may be too small for robots to reach, or the robot may damage sur¬ 
rounding tissue if it enters the region. An alternative is for the robot to extend a narrow probe into 
the region, e.g., to measure chemicals or release drugs. For instance, sensors, tens of nanometers 
in size, can be inserted through cell membranes 10,15 
micron-scale robots considered here. 


These probes are much smaller than the 


Another application for small probes arises in biological tissues with size-dependent viscosity 28 


In such tissues, micron-size objects experience larger viscosities than smaller ones. In such cases, 
nano-scale probes could extend through regions of low viscosity that are too small to accommodate 
the whole robot. 


4.1 Geometry 

To focus on the motion of the probe itself, we assume the robot is anchored to a larger structure, 
so the robot does not move as the probe extends. For simplicity, we consider a probe extending 
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Figure 3: Telescoping probe with n = 5 segments, (a) Probe extending from a spherical robot, (b) 
Cross section of the probe with an exaggerated vertical scale. 










h 


telescoping probe 



Figure 4: Coefficients for viscous drag, /ifl u id> an d internal dissipation, /linternab vs. extent of tele¬ 
scoping probe with n = 5 segments, r/L = 0.05 and s/L = 0.02. 


perpendicular to the robot surface, giving the axisymmetric geometry shown in Fig. [3] 

With / giving the fraction of each segment’s extension, we consider uniform extension so / = 
F/T , and the probe extends from / = 0 to / = F in time T. Convenient choices for characteristic 
size and speed are the segment length and speed of the probe’s end: 

d 0 = L (13) 

Vo = (n-l)Lf (14) 


4.2 Behavior 


We determine for the telescoping probe using the same simplifications as described above for 

expanding robots, i.e., consider the external surface of the probe to be truncated cone. The end 
of the cone moves at speed (n — 1 )L/, and the speed along the side interpolates linearly from 0 at 
the base to (n — 1 )Lf at its end. The result, shown in Fig. [ 4 J shows larger variation than for the 
expanding robot. This is due to the narrower width of the probe, so a larger proportion of drag 
arises from the extending sides compared to the constant area of the tip. 

Sliding surface area between segments i and i — 1 is Si = 27r(r + [n — i)s)L{ 1 — /). The surfaces 
slide past each other with relative speed Vi — i>i-i = Lf. Internal power use, from the sliding friction 
between adjacent shells, is similar to that of the expanding robot: 


-^internal — E kf 

riction*^i (Lf) 2 


(15) 


With Eq. [6j this gives 

^internal = 7r(l - f) — — ((n - 2) y + 2y) (16) 

n-1 \ L LJ 

Thus /^internal decreases as the probe expands, i.e., as / increases, since each segment is in contact 
with less of its neighboring segments. Fig. [4] shows the behavior of ^internal- 

4.3 Example 

Table [3] illustrates the situation for the two scenarios of Tabic [lj The Reynolds and Womersley 
numbers are small in both cases. Combined with the geometric coefficients shown in Fig. [4j these 
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scenario 


LOW 

HIGH 

time to extend 

T 

1 ms 

Is 

speed of extending tip 

vo 

2 mm/s 

2 /im/s 

Reynolds number 

Re 

2 x 10" 3 

2 x 10~ 10 

Womersley number 

W 

0.03 

10" 5 

fluid power factor 

rjdoVo 

4 x 10 -3 pW 

4 x 10“ 5 pW 

internal power factor 

^friction^o^O 

4 x 10" 3 pW 

4 x 10" 9 pW 

energy used 

E 

6.8 x 10" 18 J 

6.4 x 10“ 17 J 

glucose molecules 


3 

30 

Table 3: Telescoping probe for the scenarios of Table [l] with do = L = 

1 /im, r = 50 nm, s = 20 nm, 

n = 5 segments, and fcf r i c tion = 10 3 kg/ ( 

m 2 s) for a probe extending at 

a uniform speed from / = 0 

to f = F = 0.5. The energy listed is for expanding at 

a uniform rate, 

given both in joules and the 


number of glucose molecule oxidations required to produce that energy, at 50% efficiency 20 
both scenarios, fluid viscosity accounts for over 90% of the energy used. 


In 


values indicate viscous drag is the dominant contribution to overall power use, and the power required 
is well below 1 pW. 

Minimizing energy requires extending the probe about twice as fast at the beginning of the 
extension than at the end. This only reduces energy by about 2% compared to uniform expansion. 
Of potentially greater significance is balancing power use within the robot: the energy minimizing 
extension uses power uniformly, at a rate 40% lower than the maximum required during uniform 
expansion (which occurs as the probe reaches its maximum extent). Comparing with the smaller 
improvement for the expanding robot shows that a nonuniform rate of shape change is mainly 
beneficial for extending long, narrow structures, where the total area exposed to fluid drag changes 
significantly during the shape change. 


5 Aggregating Robots 


In the above examples, a single robot changes its shape. Another type of metamorphic robot is a 
modular robot, which is an aggregate of smaller robots (i.e., modules) that move with respect to 
each other. The aggregate changes shape primarily due to motions of its constituent robots rather 
than shape changes by those constituents, although they may adjust their shapes to fit into the 
aggregate. This section focuses on the behavior of the smaller robots, i.e., the modules, and takes 
the aggregate to be a passive structure. 


In some cases, modules remain connected in an aggregate as they move 23 43 


In other cases, 

more relevant for robots in fluids, the modules can move apart and congregate as needed. Such 
swarms of separate modules can form aggregates with variable numbers of modules, even with simple 

00 


controllers and noisy sensors 


This involves navigation for the modules to find each other, e.g., 


using chemical signals as used by microorganisms (3 12 or acoustic signals 21 37 . Hydrodynamic 
interactions could also aid coordination of nearby robots 33 . This type of modular robot introduces 


a new ingredient not found in the previous examples: robots moving through the fluid toward each 
other. 


For simplicity, we consider robot motion in an otherwise stationary fluid and take the robots 
to be neutrally buoyant spheres. This discussion does not include dissipation associated with any 
locking mechanisms when modules attach to the growing aggregate, or arising from small shape 
adjustments, if any, the modules require for this attachment. 
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5.1 Geometry 


To form an aggregate, the robots propel themselves toward each other through the fluid. Power 
dissipation depends on their locomotion method. Methods involving extended structures, such as 
flagella, pose a problem of tangling as the robots approach each other. To avoid this problem, we 
suppose robots use tangential surface motion for propulsion 130 . Specifically, we consider treadmills 
in an equatorial band covering half the robot surface area (T9j. With this method, the robot does 
not change shape as it moves, so internal power dissipation depends on the speed of the treads but 
not on the position of the robot relative to other robots. 

In the previous examples, robot actuators directly altered the shape of the robot, characterized 
by the parameter /. In the case of modular robots, the configuration is the relative positions of the 
modules. Actuators do not directly specify the configuration. Instead, actuators exert force against 
the fluid to move the robot. The speed of this motion depends on the positions of other robots. Thus 
there is not a fixed relationship between actuator motion and configuration of the aggregate. This 
means robots moving in fluid are a more complicated shape-changing situation, with two properties 
determining the behavior, i.e., the relations of tread speed to locomotion speed, and of tread speed 
to fluid dissipation. 

Convenient choices for characteristic size and speed are the sphere radius and tread speed over 
the sphere surface: 


d 0 = a (17) 

^0 = ^tread (18) 

In addition to characterizing power dissipation (Eq. [3] and Eq. [6]), the tread speed determines 
the sphere’s locomotion speed u according to 19 

U = I^locTtread (19) 

where the proportionality h\ oc varies with the sphere’s position. 

To illustrate the behavior of modular robots, we consider two cases, shown in Fig. [5j First, a 
sphere approaches a much larger, stationary aggregate, modeled as a sphere approaching a planar 
wall. In the second case, two identical spheres approach each other along their common axis. This 
illustrates the behavior at an early stage of aggregation, i.e., when the first two robots approach 
each other to start forming an aggregate. 

We characterize the position of the sphere by the distance L of the center of the sphere to the 
wall, in the first case, and from the center of one sphere to the midpoint between the two spheres, 
in the second case, so the distance between their centers is 2 L. We define the normalized distance 

6 = L/a (20) 

so 5 = 1 corresponds to the sphere touching the wall or the other sphere in the two cases, respectively. 

5.2 Behavior 

We determine h\ oc and ha u ^ for the sphere by assuming quasi-static Stokes flow. 

Fig. [6] shows /Hoc for the two cases. As the separation decreases toward zero (i.e., S —» 1), the 
treadmill locomotion becomes increasingly less effective, i.e., a given tread speed produces slower 
motion through the fluid. With two approaching spheres, the distance between them decreases at 
twice the speed of each sphere with respect to the fluid. 

We consider a situation where the robots start relatively far apart, with 5 = 5q somewhat larger 
than one, and move toward each other to a smaller distance <5i ss 1. Since h\ oc —> 0 as L —> a, we 
consider motion of robots to within a small positive distance, e.g., 50 nm, which is sufficient for the 
robots to link to each other. This corresponds to Si = 1.05 for a 1 /rm-radius sphere approaching a 
wall, and S\ = 1.025 for two spheres approaching each other. 
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Figure 5: Cross section of a sphere approaching a wall or another sphere. 



Figure 6: Ratio of locomotion speed to tread speed, h \ oc , vs. normalized distance 
for a sphere approaching a wall and for two spheres approaching each other. 


8 = L/a (Eq. 20) 


h 


sphere 



25 - 


20 - 


L / a 


Figure 7: Coefficients for viscous drag, ft fluid (solid), and internal dissipation, /internal (dashed), 
vs. normalized distance 8 = L/a (Eq. 20) for a sphere approaching a wall and for two spheres 


approaching each other. Internal dissipation is the same for both cases and independent of distance. 
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scenario 


LOW 


HIGH 


tread speed 
Reynolds number 
Womersley number 
fluid power factor 
internal power factor 
diffusion coefficient 


Re 

W 

Vdovo 

^friction 

ksd body/(b'7n7^sphere) 


lmm/s lO^tm/s 

10-s 10 -g 

0.03 3 x 10 -5 

io- 3 P w io- 3 P w 

io- 3 P w io- 7 P w 

10 " 13 m 2 /s 10 -17 m 2 /s 


Table 4: Sphere motion for the scenarios of Table a with do = a S phere = lM m and ^friction = 
10 1 * 3 kg/(m 2 s) for a sphere with uniform tread speed in an equatorial band covering half the surface 
area, ks is the Boltzmann constant and Tbody = 310 K is body temperature. 


scenario LOW high 


case 


wall 

two spheres 

wall 

two spheres 

final normalized distance 

Si 

1.05 

1.025 

1.05 

1.025 

motion time 

T 

10 ms 

9 ms 

Is 

0.9 s 

energy used 

E 

4.9 x 10" 16 J 

4.6 x 10" 16 J 

3.0 x KT 14 J 

2.8 x 10" 14 J 

glucose molecules 


250 

230 

15000 

14000 


Table 5: Sphere motion for the scenarios of Table |4j starting from normalized distance <5o = 5 with 
uniform tread speed. Each scenario shows behavior for two cases: a sphere approaching a wall and 
two spheres approaching each other. The energy listed is given both in joules and the number of 
glucose molecule oxidations required to produce that energy, at 50% efficiency [20 . In the LOW and 
high scenarios, fluid viscosity accounts for 60% and nearly 100% of the energy used, respectively. 


Fig- 0 shows the power dissipation: values for ha u id and hi nte rnai for the two cases. The internal 
dissipation coefficient is independent of distance. When the distance is a few multiples of the sphere 
radius, fluid power dissipation is close to that of an isolated sphere using this tangential motion^ 


5.3 Example 


Table |4| illustrates the behavior, including diffusion coefficients of spheres [2], for the scenarios of 
Table[lJ Table[5]gives the time and energy used for the two scenarios. Over these times, the diffusion 
coefficients given in Table [4] indicate the effect of Brownian motion: moving the spheres by typical 
distances \f&DT of about 100 nm and 10 nm, for the low and high scenarios, respectively. These 
distances are small compared to the range of motion considered here, i.e., several microns. Thus, 
Brownian motion is noticeable but fairly minor for the scenarios considered here. 

The tread speeds in Tabic [4] and geometric coefficients in Fig. [7] indicate viscous drag is the 
dominant contribution to overall power dissipation, which is well below 1 pW for the times in Table[5] 
The dissipation is per sphere, so total power in two-sphere case is twice the value given here. The 
relatively small variation in dissipation as a function of position (Fig. [7| leads to only a slight 
reduction in energy use (0.1% or less) when optimizing the tread speed as the robot moves. 


An additional issue for aggregating robots arises if they obtain power from chemicals in their 
environment. In that case, as robots approach there is less power available due to competition for 
those chemicals from other robots 31 . This could increase the time required to aggregate many 


robots in small volumes. However, this is not a significant issue for modest numbers of robots 
obtaining power from oxygen and glucose dissolved in the fluid. This is because, for the scenarios 
in Table [5j the energy requirements are considerably less than available chemical power for a few 
dozen robots operating in close proximity in or near blood vessels 


20 


1 The dissipation shown here is somewhat larger than the value ignoring fluid vorticity [l~9| . due to additional 

dissipation arising from the velocity gradient at the ends of the treadmill. Dissipation due to this edge-effect depends 

on the distance at the ends of the tread over which velocity changes. This study uses a distance corresponding to 

50 nm bearings for the treadmill |19|. 
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6 Discussion 


Microscopic metamorphic robots with the parameters studied here will likely have enough power 
from chemicals available in biological tissues 13 20 . Biological cells indicate this is physically 


possible, since cells actively change their shapes in spite of facing similar power limitations and 
dominance of surface forces as the robots considered in this paper. 


There are significant engineering challenges to creating these robots. These include building 
robots from materials whose sliding friction is as low as assumed with Eq. [5] With this value, 
dissipation from internal friction and external viscous drag are of comparable magnitude for the 
examples considered in this paper. Thus there is is relatively little decrease in power use from 
additional reduction in friction. On the other hand, building the robots from conventional materials, 
which have much larger friction, will lead to significantly more power dissipation, almost all of which 
will be due to internal friction. 


The focus in this paper is on shape change involving a single degree of freedom. More complex 
changes include extending multiple probes or manipulator arms. Similarly, modular robots can 
involve several modules moving simultaneously on or near the surface of the aggregate |4 24 34 


Simultaneous moving parts lead to complex hydrodynamic interactions, similar to those for nearby 


objects moving in a fluid 18 22 29 


Extensions to this study could examine other types of reconfiguration. For instance, because 
surface forces dominate behavior of microscopic robots, changing surface properties such as affinity 
for water 


27 could complement shape changes in allowing robots to adapt to their tasks. 


A final direction for study is developing controls suited to microscopic metamorphic robots, 
particularly their limited computation and noisy sensors. One approach, for modular robots, is 
extending simple local rules that grow structures 4,j36 j40 and other swarm algorithms [5] to account 
for large surface forces and Brownian motion. Automated rule generation 
these approaches. 


26 could complement 
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